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ABSTRACT 


When  a  metal  matrix  composite  is  cooled  down  to  room  temperature 
from  the  fabrication  or  annealing  temperature,  residual  stresses  are  induced  in 
the  composite  due  to  the  mismatch  of  the  thermal  expansion  coefficients 
between  the  matrix  and  fiber.  A  method  can  be  derived  for  calculating  the 
internal  stresses  appearing  in  Metal  Matrix  composites  of  A1  matrix  with  SiC 
particles  due  to  differences  in  thermal  expansion  coefficients.  Special  attention 
is  paid  to  creep  deformation  in  the  matrix  phase.  The  analysis  shows  that 
considerable  internal  stresses  and  creep  deformation  appear  in  the  composites 
when  subjected  to  cooling. 
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1.  INTRODUCTION 


A.  METAL  MATRIX  COMPOSITES! MMCs) 

Discontinuously  reinforced  Metal  Matrix  Composites(MMCs)  represent  a 
group  of  materials  that  combine  the  strength  and  hardness  of  the  reinforcing  phase 
with  the  ductility  and  toughness  of  the  matrix.  Aluminum  alloys! Al)  reinforced 
with  Silicon  Carbide!  SiC)  in  particulate,  platelet,  or  whisker  form  and  fabricated  by 
powder  metallurgy  methods  are  receiving  a  great  deal  of  attention  from  researchers 
and  engineers. 

B.  THE  ANALYSIS  OF  MMCs(  Al/SICs)  IN  THE  STABLE  MEMBER 

Previous  research  has  evaluated  the  dimensional  stability,  and  the  thermal 

and  mechanical  properties  of  several  Al/SiC  MMCs  in  stable-member  applications 
for  missile  inertial  guidance  system.  The  results  reveal  that,  although  the  candidate 
materials  from  a  powder  blend  of  SiC  and  A1  alloy  consolidated  by  Vacuum  Hot 
Pressing! VHP )  into  cylindrical  billets,  followed  by  Hot  Isostatic  Pressing! HIP)  for 
full  densification,  have  better  isotropic  properties,  these  MMCs  show  some  micro 
creeping  behavior  in  service.  The  microcreep  of  these  MMCs  will  affect  the 
dimensional  stability  of  stable  members,  and  is  likely  to  arise  from  two  sources:  ( 1 ) 
creep  of  the  metal  matrix  caused  by  internal  stresses  creep  conditioned  by  externally 
applied  stresses  and/or  (2)  phase  transformations  during  the  creep  condition.  A  first 
step  toward  understanding  the  cause  of  the  dimensional  stability  problem  is  to 
analyze  the  influence  of  internal  stresses.  The  internal  stresses  are  thermal  residual 
stresses  when  a  Metal  Matrix  Composite  is  cooled  down  to  room  temperature  from 


the  fabrication  or  annealing  temperature.  These  residual  stresses  are  induced  in  the 


composite  due  to  the  mismatch  of  the  thermal  expansion  coefficients  between  the 
aluminum  alloy  matrix  and  the  silicon  carbide  particles.  The  model,  based  on 
Eshelby's  for  mismatch  problems  with  simplified  linear  elastic  material  behavior  for 
both  particles  and  matrix,  has  been  used  to  solve  the  problem  of  thermal  residual 
stress  [ref  1]. 

C.  PURPOSE 

The  purpose  of  this  thesis  is  to  calculate  the  thermal  residual  stress  of  Al/SiU 
using  theoretical  methods  and  then  to  determine  what  influence  this  thermal 
residual  stress  has  on  creep  deformation.  In  order  to  do  this,  the  thermal  residual 
stresses  will  be  determined  by  focusing  on  elastoplastic  matrix  and  elastic  inclusions 
rather  than  on  elastic  matrix  and  elastic  inclusions.  Therefore,  the  case  where  both 
phasesf  matrix  and  inclusions  )  are  perfectly  elastic  will  be  treated  first.  Then,  the 
theoretical  results  will  be  compared  with  previously  obtained  experimental  data  for 
Al/SiO  composites. 


II.  LITERATURE  REVIEW 


A.  MMCs 

Metal  matrix  composites! MMCs).  including  eutectic  composites,  are  becoming 
important  in  applications  as  structural  components  which  are  to  be  at  intermediate 
and  high  temperatures.  When  MMCs  are  fabricated  at  high  temperature  or  annealed 
at  certain  high  temperatures,  the  MMCs  have  undesirable  properties,  such  as  low 
tensile  yield  and  ultimate  strengths.  Those  results  are  mainly  due  to  residual 
stresses  that  are  caused  by  the  mismatch  of  the  thermal  expansion  coefficient 
between  the  matrix  and  fiber.  The  residual  stresses  so  induced  have  been  observed 
in  tungsten  fiber/copper  composites  and  in  SiC  whisker/6061  A1  composites  and 
were  based  on  analyzed  a  one-dimensional  model  for  a  continuous  fiber  system  or 
spherical  particle  system.  The  model,  based  on  Eshelby's  equivalent  inclusion 
method,  has  been  used  to  solve  the  problem  of  thermal  residual  stresses[ref 

B.  ELASTIC  MATRIX  AND  ELASTIC  INCLUSIONS 

Let's  consider  a  two— phase  suspension  of  finite  extent  in  which  one  phase  is  a 

matrix  while  the  other  is  in  the  form  of  inclusions.  It  is  assumed  that  the  inclusions 

can  be  treated  as  identical  ellipsoids  with  corresponding  axes  aligned.  Further,  the 

inclusions  are  assumed  to  be  randomly  distributed  in  the  matrix  in  such  a  way  that 

the  suspension  is  homogeneous  on  a  macroscopic  scale.  The  task  of  the  present 

section  is  to  solve  Eshelby's  transformation  problem  for  the  case  where  the 

fractional  volume  f  of  the  inclusions  is  finite.  Following  Eshelby  [ref  2],  we  consider 

first  an  auxiliary  problem,  what  is  the  resulting  stress  a  ^  in  an  inclusion  when  every 

ij 

inclusion,  having  isotropic  plastic  moduli  equal  to  those  of  the  matrix,  is  subjected 


T 

♦  o  a  uniform  transformation  strain  6  17  To  find  the  exact  solution  appears  to  he 

*  l 


impossible  since  no  detailed  information  is  given  about  the  spatial  distributi 


ion  ol 


the  inclusions.  All  we  may  find  in  this  situation  -'re  approximate  average  fields.  Tlu 

average  constrained  strain  €  ~  in  anv  inclusion  may  be  written  in  the  form 

i  j 


i  n 

6-C  =  6C  +  6C 


1) 


'  J 


‘  J 


(  1  ) 


C 

Here,  6  is  the  constrained  strain  in  a  typical  inclusion  which  would  be  produced  if 
>  i 

it  alone  is  embedded  in  an  infinite  elastic  body.  This  term  is  given,  according  to 
Eshelby[ref  ].  as 


€C=S  €T 

l]  ijkl  i) 


(  2  ) 


Where  S  is  Eshelbv's  tensor  which  depends  on  the  aspect  ratios  of  the  ellipsoidal 

i]  kl 

inclusion  and  Poisson's  ratio. (Here  and  in  the  following  we  use  the  usual  summation 

1 1 

convention:  the  range  of  subscripts  is  form  1  to  3).  is  a  contribution  due  to  all 

* ) 

the  surrounding  inclusions  and  the  presence  of  free  boundary;  this  term  is  entirely 

elastic,  and  is  assumed  to  be  effectively  constant  in  and  around  the  typical 

1 1 

p 

inclusion.  Then,  it  will  be  a  good  approximation  to  regard  €*  as  the  average 

»  i 

elastic  strain  in  the  matrix  for  the  contribution  of  any  one  inclusion  to  the  average 
field;  outside  it  is  negligible  when  there  are  a  great  number  of  inclusion  in  the 
material.  For  a  homogeneous  elastic  body  in  static  equilibrium,  a  volume  integral  of 
any  component  of  the  elastic  strain  associated  with  an  internal  stress  should  vanish 
provided  that  the  integration  is  taken  throughout  the  body.  Thus 


V.  FJT  HTUUl  WWJ 


j.i  ijAt  »ii  .■*  >.i  l.>  ».i^«-i*i.>  ».i  wL«tl.»d.«i 


C  -C  -T 

1-f)  €  -  +  f ( €  '  —  €  1  )  =  U 

*  J  ‘ )  >  J 


From  ( 1),  (2)  and  ( 3),  we  have 


€4:  =  (l-ns  €  J  +  f€T 

i  j  ijkl  k  I  l  ) 


The  average  interna)  stress  in  any  inclusion  is  given  by 


P  =  C  (€-C-£T) 

ij  ijkl  ki  kl 


Where  C  is  the  elastic  moduli  ol  the  matrix  (and  also  of  the  inclusions  in  the 
ijkl 

present  situation)  [ref  3].  Where  the  inclusions  have  elastic  moduli  differing  from 

those  of  the  matrix,  the  same  approach  as  used  by  Eshelby  [ref  2]  will  be  applicable. 

* 

The  essence  of  this  approach  is  terms  of  equivalent  inclusions.  Let  C  be  the 

i  j  kl 

T* 

elastic  moduli  of  the  inclusions,  and  €  be  a  uniform  transformation  strain  of  everv 

'].  T 

inclusion.  Then  the  transformation  strain  of  the  equivalent  inclusions,  6  ,  can  be 

1 J 

determined  by  solving  a  set  of  equations 


P  T  *  C  T* 

C  (  ’  -  e  1  )  =  C  (  €-  -  -e1  ) 

ij  kl  k  1  k  1  ijkl  k  1  k  1 


with  (  4  ).  The  average  stress  in  any  inclusion  is  given  by  (  5  ),  of  course.  Consider 
a  composite  specimen  which  is  free  of  stress  at  a  references  temperature  To-  Let  the 

matrix  and  the  inclusions  have  isotropic  linear  thermal  expansion  coefficients  a  and 

* 

a  ,  respectively.  The  elastic  state  when  the  specimen  is  subjected  to  a  uniform 


& 


¥ 


temperature  change,  aTq=T— To,  is  described  by  the  situation  where  the  matrix 


and  the  inclusions  undergo  uniform  transformation  strains  oaT^  and 

♦  ,  11 
a  aT£  .respectively.  Here  6  is  the  Kronecker  delta.  The  resulting  state  of  stress 

‘j  '  *) 

can  be  analyzed  by  putting 

€T*  =  (a*  —  a)  aT  6  (7) 

>  j  »i 

The  overall  thermal  expansion  coefficients  a  can  be  expressed  as  a  function  of  a. 

*  * 

a  ,  C  ,  C  ,  f  and  the  aspect  ratios  of  the  inclusions. 

ijkl  ijkl  r 

C.  ELASTOPLASTIC  MATRIX  AND  ELASTIC  INCLUSIONS 

When  the  temperature  change  is  great,  one  should  expect  plastic  deformation, 
either  in  the  matrix  and/or  in  the  inclusions,  which  relaxes  the  internal  stress 
accumulated.  An  important  practical  case,  the  case  where  the  plastic  deformation 
occurs  only  in  the  matrix,  will  be  considered  in  the  following.  Further,  for 
simplicity,  we  assume  that  the  plastic  strain  is  uniform  throughout  the  matrix  and 
the  plastically  nondeforming  inclusions  can  be  treated  as  identical  spheroids.  Then, 
the  plastic  strain  in  the  matrix,  6  P,  may  be  written  in  consideration  of  a  symmetric 
mode  of  plastic  deformation  and  the  volume  constant  law,  as: 

€p  =  €  \6  6  -hs  6+6  6)]  (  8  ) 

i)  P  3i  3j  ^  li  Ij  2i  2) 

** 

Where  €D  is  a  temperature— dependent  parameter.  In  this  case,  the  overall  strain  € 

—  T 

of  the  composite  is  given,  with  use  of  Equation  £ .  =f€  ,  by 


p 

$ 

! 


*:* 

I*! 

;% 

u5 


a 

Sj 

a 

Sf 

3 


$ 

1 

*1*1 


a 


€  =  oaT f>  +€  P+fG  1 

*J  *]  >J  >1 


Where  €  J  is  calculated  from  (4)  and  (6)  with 
>  j 


=  (a  —  a)AT<$  — €  P 
>j  ‘i  >i 


Now,  the  temperature  dependence  of  £p  will  be  found  through  an  energy— balance 
consideration.  Since  energy  dissipative  process  is  involved  in  plastic  deformation, 
the  decrease  in  the  elastic  energy.  —  <5Eel ,  corresponding  to  a  virtual  plastic 
deformation  6€p,  should  be  balanced  against  the  dissipation  energy.  If  the  matrix 
material  has  a  constant  flow  stress  f  in  a  uniaxial  deformation,  the  condition  is 

y 

written  as: 


Where  Eei  =  — 

=  Ci€p  +  (-2^p(  a  —  o)aT 
+  CM  a*  -  oM2 


v3 


It  is  noted  that  when  the  composite  is  elastoplatic  the  following  condition  should  be 


satisfied: 


(Eel)j  >  (Eel)jj 


£ 

,"A 


m 

Wi 


Where  the  subscripts  I  and  II  refer  to  the  entirely  elastic  case  and  the  elasfoplastic 
case,  respectivelyfref  3]. 

D.  CREEP 

If  a  stress  is  suddenly  applied  to  pure  metals,  some  solid  solutions  and  most 
engineering  alloys  at  a  temperature  near  or  greater  than  0.5Tm( where  Tm  is  the 
absolute  melting  point  of  the  metal  or  alloy)  deformation  proceeds  as  shown  in  Fig. 
2.1.  The  initial  application  of  stress  causes  an  instantaneous  elastic  strain  €o  to 
occur.  If  the  stress  is  sufficiently  high  ar>  initial  plastic  deformation  6P  also  occurs. 
At  low  temperatures,  significant  deformation  ceases  after  the  initial  application  of 
stress  and  an  increase  in  stress  is  required  to  cause  further  deformation.  At  elevated 
temperatures,  deformation  under  a  constant  applied  load  continues  with  time.  The 
early  stage  of  such  deformation  called  primary  creep  is  characterized  by  an  initially 
high  creep  rate  de/dt  which  gradually  decreases  with  time.  Eventually  a  linear 
variation  of  creep  strain  accumulation  with  time  is  observed.  This  steady— state 
creep  region  is  characterized  by  a  constant  minimum  creep  rate.  Steady— state  creep 
rates  depend  significantly  on  stress  and  temperature  and  are  used  frequently  to 
compare  the  creep  resistance  among  alloys.  After  significant  deformation  in 
steady— state  creep,  necking  occurs  or  sufficient  internal  damage  in  the  form  of  voids 
or  cavities  accumulate  to  reduce  the  cross-sectional  area  resulting  in  an  increase  in 
stress  and  creep  rate.  The  process  accelerates  rapidly  and  failure  occurs.  This  region 
of  the  creep  curve  is  called  tertiary  creep.  Fig.  2.1  also  shows  the  derivative  of  the 
creep  curve  or  the  creep  rate  curve.  Fig.  2.1  represents  the  most  common  creep 
strain/time  behavior  and  Fig.  2.2  illustrates  the  strain  rate/time  behavior  at 
temperatures  below  approximately  0.3Tm.  In  Fig.  2,  only  transient  creep  is 


observed  which  is  characterized  by  a  continuously  decreasing  creep  rate  that 
approaches  zero  as  the  inverse  c>f  time.  Such  low— temperature  crpep  behavior  is 
called  logarithmic  creep.  Fig.  2.2  also  illustrates  the  creep  rate/time  behavior  of 
alloys  that  exhibit  a  continuously  increasing  creep  rate  in  the  early  stages  of  creep. 
Such  alloys  have  been  designated  Class  1  alloys.  Creep  by  viscous  dislocation  glide 
results  when  dislocations  glide  or  move  in  slip  planes  under  the  action  of  an  applied 
stress.  These  dislocations  drag  along  solute  atoms  attracted  to  the  strain  fields  of 
the  dislocations.  In  order  for  the  dislocations  to  move  the  solute  atmospheres  must 
diffuse  in  the  direction  of  dislocation  motion.  At  lower  temperaturesf  ~0.3Tm)  or  at 
higher  stress  levels,  creep  occurs  by  thermally  activated  dislocation  glide.  Under 
these  conditions  dislocations  can  overcome  barriers  to  motion  without  dislocation 
climb.  The  localized  motion  of  short  segments  of  small  dislocations  is  important  in 
overcoming  barriers.  Because  of  the  larger  contribution  of  stress  in  thermally 
activated  dislocation  glide  a  temperature  and  stress  dependence  for  creep  different 
from  that  for  dislocation  climb  or  viscous  glide  will  be  observed. fref  61 


Creep  rarte.  dc/dt  Creep  strafn.  e, 
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Fig.  2.2  Creep  rate  vs  time 
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Fig.  2.3  Typical  creep  curve 
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III.  THEORETICAL  BACKGROUND  AND  DERIVATION 


A.  THEORETICAL  BACKGROUND 

To  simulate  Metal  Matrix  Composites!  MMCs),  we  consider  an  infinite 
isotropic  body  containing  an  ellipsoidal  fiber  as  shown  in  Fig.  3.1.  The  infinite 
body! composite)  experiences  a  temperature  T1  —  TO.  where  T1  is  assumed  to  be 
larger  than  the  relaxation  temperature  of  the  matrix  Tr  and  of  course,  Tr>T0.  To 
make  the  problem  simpler,  the  case  of  one  fiber  embedded  in  the  infinite  body  is 
considered  here. 


Fig.  3.1.  Theoretical  model 

Based  on  the  model  in  Fig  2,  let  us  consider  ♦*»<»  rase  **  elastic 
inclusions— elastic  matrix,  elastoplastic  matrix— elastic  inclusions  and  creep 


11 


deformation: 


1.  Elastic  Inclusions  and  Matrix 

-C 

First  let  us  consider  the  average  internal  constrain  strain  €  .  To  find 

_C  .  ,J 

6  ,  we  insert  Equation(l)  and  Equationf-3)  into  Equation(4)  and  find 


e~^=(  1— f)S  e  T  +  f€  T, 


€  =(1— f)S  G  +  f€  ,  which  represents  the  average  internal  constrain  in  anv 

i  j  ij  kl  k  1  i  j 

inclusion  in  terms  of  transformation  strain.  Next,  the  average  internal  stress  in  any 

T 

inclusion  can  be  represented  by  Equation(5).  Here,  6  is  an  unknown  term. 

» j 

Therefore,  to  find  this  unknown  term  we  must  consider  Equation  (6)  and  Equation 

T* 

(1).  In  Equation  (6),  *  represents  the  inclusion,  and  €  represents  the  uniform 

11  .  T 

transformation  strain  for  every  inclusion.  Once  we  have  determined  €  ,  according 

T  1 J 
to  the  Equation  6  =  fg  1  we  can  determine  the  overall  strain  of  the  specimen,  and 
ij  ij 

T* 

in  accordance  with  Equation! 7)  we  can  determine  the  uniform  transformation  6 

i  j 


2.  Elastoplastic  Matrix  and  Elastic  Inclusions 

The  case  for  elastoplastic  matrix— elastic  inclusions  is  different  from  that 

of  elastic  inclusion  and  matrix.  This  is  due  to  the  fact  that  in  this  case  plastic 

deformation  occurs  only  inside  the  matrix,  so  we  must  consider  the  plastic  strain  in 

the  matrix.  Therefore  we  can  find  €  ^  in  accordance  with  Equation(o)  and 

ii 

1  I  i  ^  I 

Equation!  10).  From  Eej  =  —  ^f<x  t  we  can  represent  a  =F(€  and 

-  ij  ij  ij  ij 

T*  n 

6  .  =F(€.f).  After  representing  them  this  way, 


v  h  -IeT* 

Eel  =  -  nf<r  6 

t  2  2 

=C1  €p2  +  C2  €  ( a*  -  a)AT  +  C2  €  ( a  -a)  *T  (  14  ). 
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Differentiating  this  with  respect  to  €  ,  we  can  obtain  the  value  of  as  shown 


6  =  -o^y  ( a  ~  a)&!T  +  f)<r 


(  15  ) 


Therefore,  using  the  value  of  €  to  find  the  value  €  the  value  of  cr'^  can  be 

P  ij  ij 


obtained,  and  inserting  this  value  into  Equation(16),  we  are  ultimately  able  tc 
obtain  th°  total  average  stress  in  the  matrix. [ref  9] 


<  c r-1  >  =  C-M  1  <  €  >,-  =  -f 
ij  ra  ijkl  ki  VR  ij 


3.  Creei 


Applying  the  value  of  the  internal  stress(above)  to  Equation(13).  we  can 


determine  creep  deformation.  Following  from  this  theoretical  base,  let  us  first 
calculate  mathematically  the  thermal  residual  stress  of  Al/SiC  MMCs. 


B.  DERIVATION 


1.  Elastic  Inclusions  and  Elastic  Matrix 


First  let  us  calculate  €  1  for  the  case  of  elastic  inclusions  and  elastic  matrix. 

T  ,J 

In  order  to  determine  €  1 ,  we  rewrite  Equation(4)  in  matrix  form: 

i  i 


S  1 1  S  (2  S  1 3  S  1 4  S  15  S 1 6 
So  I  S22  S23  S24  S25  S26 
-n_n  S31  S32  S33  S34  S35  S36 
S,,  S42  S43  S44  S45  S46 
S51  S52  S53  S54  S55  S56 
.  Sgl  S 62  S63  S 6 4  S 65  See 


+  t  !s 

€44 


it.  ika  li.l'i  t’i  d'i  I'k, 


PM 
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Simplifying  this  matrix  we  obtain: 


fit  S 1 1  f  +  S  i  2  €.i  +  S1363I 
S2 1  €  j  ^  +€o^(  So  2  +  f  )  +  S  2  3^3^ 


— n  ff  j €  1 1  +  S 3 2 €2^  +  £3^  S33  +  f  )  ,  /  €33 

— v1— c  t/ o.  ,  ,  x  \  +  1  c.t 


€ 4 4 (  S44  +  f  ) 

€5§( S55  +  f  ) 

€66(  Se6  +  f  ) 


Further,  simplifying  this  matrix  we  find: 


£m  =  GiK  1 — f)Su  +  €2^12(1— f)  +  63$Si3(l— f)  4-  f€,i 
€22  =  €ii(l-f)S2!  +  €2^22(1— f)  +  £3^23(1—1)  +  f£2$ 
£33  =  £iH  1 — f)Ssi  +  £2lS32(l— f)  +  £3^33(1—0  +  {£3$ 
£44  =  €4$(l-f)S44  +  f£4$ 

€55  =  €sS(l— f)Ss5  +  {£5$ 

£66  =  £6^1“ : 0$66  +  f£6^ 


(  17) 


Refer  to  Appendix  (A)  for  the  values  of  Sijkl- 
Next,  let  us  rewrite  Equation(6)  in  matrix  form. 


(l-i/)  v  i/OOO 
v  (l-i/)  v  0  0  0 

v  v  ( 1-v )  0  0  0 

0  0  0  2n  0  0 

0  0  0  0  2/x  0 

0  0  0  0  0  2/i 


(l-i/)  v  j/000 
v  (1-*/)  v  0  0  0 

v  v  (1-*/)  0  0  0 

0  0  0  2p  0  0 

0  0  0  0  2/i  0 

0  0  0  0  0  2/i 


Simplifying  the  above  matrix,  we  have  the  following; 
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€u(l— +  I'Geji+J'Gcsjj— (1— t^Hei) 

—  MeJ  —  ^Hesi 

*/)J+€js^J=i<3€^  +  (1— t/)G€2i+t/G€s&— — 
(1—  y)H€2J  —  ^HGsJ 

€ni/J+€jjX'J+€ss{l— i/)J=x>G£i)+  jXj£2J+(1— i/)G€ji— j/H€i!  — 


*He2$  -  (i-^)H€SS 


=  o 

€55  =  0 
^66  =  0 


From  Equation  (17)  and  (18),  if  we  simplify  and  solve  for  6  ,  we  obtain: 

» j 


K€|)+L62i+M€3^='— (1 — 
Nei)+Pe25+Qe3S=-i/Hei)  — (1—  i/)Hc2§  — j'HGjS 
Rei)+S€2i+T€3^=— i/H€i)  — i/H€2$  — (1— 


e«$  =  e<\ 

€5$  =  €s£ 

€6^  =  tefc 


(  19) 


Now,  we  have  6  unknown  values  and  6  Equations.  Therefore,  applying  Equation(19) 


to  determine  €  : 


-YL  -  XM  -  HAVl+j/) 


m 


Putting  the  above  equations  into  simplified  form: 


i 


r  z,  1 

€2? 

Y. 

€33 

€4} 

,e6|. 

L€«| 

Here,  refer  to  Appendix(C)  for  constants  B  through  Z.  Thus,  inserting  the  value 
T  . 

€  into  Equation(5),  which  is  the  average  internal  stress  in  any  inclusion,  we  can 
determine  the  internal  stress. 


.  1 1 

tfiij 

<*22  j 
^”1 

N-* 

1 

<5 

O 

O 

O 

BZ'+CY'+DX1— Z'  ‘ 

=  G 

1/  ( l-i/)  v  0  0  0 
v  v  (l-i/)  0  0  0 

CZ'+BY'+DX'-Y' 

FZ'+FY'+EX’-X' 

<74  4  j 

<^65j 

0  0  0  2/4  0  0 

0  0  0  0  2/4  0 

0  0  0  0  0  2/4 

0 

0 

L  0  J 

2.  Elastoplastic  Matrix  and  Elastic  Inclusions 

Up  to  this  point  we  have  found  the  internal  stress  in  the  case  of  an 
elastic  matrix  and  elastic  inclusions.  Now  we  will  consider  the  case  of  elastoplastic 
matrix  and  elastic  inclusions.  Here  also  the  method  of  determining  the  internal 


and  elastic  inclusions.  Here  also  the  method  of  determining  the  internal  stress  is 

similar  to  the  case  of  elastic  matrix  and  inclusion.  However  in  this  case,  since  6  p, 

T*  'J 

the  plastic  strain,  is  included  in  the  term  6  ,  the  next  case  is  more  complicated. 

» i 

Therefore,  before  we  consider  the  relationship  between  €  p  and  6  1  ,  let  us  first  look 

T*  ij  ,J 
at  €  1  : 

cT  =(o  —  <j)aT4— €  [6  S  —  i  (6  6+5  6)]  (6  :  i= j-*l ,  i^j-*0  ).  Rewriting  the 
ij  ij  P  3i  3j  *  li  Ij  2i  2j  ij 

above  equation  in  matrix  form  we  have 


€|  *  4 1 

€22*  622 

ft  =(“*-o)AT  t 

45 

Ul  $66  . 


But  according  to  Equation(8),  €  p  can  be  shown  as  follows: 


,p  =e 


4i  4i-l/2(  4i  4i+4i  4i 
42 ~l/2f  4 2  <*12  +  42  hi 
43  4?— 1/ %  43  43+43  43 

p  4«  4«— 1/4  4  4  44+4#  44 
45  45— 1/ a  4s  45+ 4s  45 
.  4e  46— 1/2(  4s  46+46  46 


From  the  above  equation  we  can  obtain  the  following: 


€  = 


-1/2€p 

— l/2€p 

0 

0 


Therefore,  inserting  Equation(22)  into  Equation(21)  and  simplifying  we  obtain: 


9 

w 

w 

& 


'i 

'A 


i»*5 

I 

i 


V I |.».M  l.|  1.1  k|  «  | IS  l‘i  |*ft  l>»"l>i'i»i‘i*ft>ll.*>>»>WL*>4*>l 1  M  ft»i  a*4  4J  i 


T*  * 

6*  =(a  -a)AT  +  €p/2 

ej*=(a*  -  a)AT  +  Gp/2 

T*  * 

£  =(a  -  a)lT  -  Ep 

£  T*=0 
44 

£T*=0 

55 

£T*=0 

66 


Next,  using  Equation(4)  and  Equation(6),  if  6  1  is  expressed  as  €n,  the  result  is  the 

V  .  T 

same  as  Equation(17).  Inserting  Equation(17)  into  Equation(6)  and  expressing  6 

T*  *  * 

as  €  we  obtain: 

»i 


O'  rp*  m*  m* 

£n=Q|€||  —  Q2^22  —  Q3633 

T1  £p*  £p*  /Tl* 

€  =P  4€|  1  +  P562  2  +  Pe€33 
2  2 

T1  q1*  o'* 

C  =P ‘€|  '  +  P2^22  +  P3^33 

33 

£T=0 

44 

£T=0 

55 

£T=0 

66 


c  t  . 

Next,  in  order  to  find  €n,  we  insert  the  value  of  6  into  Equation(14)  and  simplify: 


rp*  m*  m* 

6u=Ri€i  |  +  H2€2  2  d-Rs^ss 

r  q1*  T*  T* 
€22=Ft4€|  1  +R$€2  2  +  R^S  3 

€«=R7€lTV  Rf^TV  R9€sT* 
€^=0 
€55=0 


Therefore  can  be  found  in  accordance  with  Equation(5).  (The  value  of  this  ax\  is 


-1 


-I 


different  from  the  value  of  <ri\  for  elastic  inclusion  and  elastic  matrix,  the  reason 


being  that  <ri|  for  elastic  inclusion  and  elastic  matrix  does  not  contain  €.P).  Here,  in 


*) 


order  to  determine  the  value  of  aj  we  must  first  find  the  value  of  €p.  From 


Gquation(lO),  if  we  differentiate  Ee]  with  respect  to  €p  we  can  obtain  Equation(ll). 

* 


Therefore,  from  ffiei  =[  2Ci€p  +  Cj(a  —  a)AT]5£p  and  —  ®ei  =  (1— f |<5€p|,  we 


can  obtain  the  value  of  6  .  Again,  inserting  this  value  of  €p  into  Equation(22)  we 


can  find  the  value  of  €.?•  Finally,  by  inserting  the  value  of  €.?  into  equation(27)  we 


»J 


obtain  the  value  of  0^.  Utilizing  Equation(lO)  again  to  find  the  value  of  € 


Eel=  —  nf 


-  1 1 
<^i'l 
<^22j 

<^33j 

*4*1 

<*55] 

^66 


€1 

^2$ 


* 


€6i 


=-Af 


X|A' 

XjA' 

XSA' 


X2€p 

X«ep 

XgGp 


A'  +  €p/2 
A'  +  6p/2 
A'  —  ep 
0 
0 
0 


=  -  \  [<#  -  X6  )  +  €PA'(  ^i+Xj+^+X^Xs+X^  ) 

+  A,2(  X,+X3+X5 )] 


(26) 


Therefore,  differentiating  Eei  in  respect  to  6p: 


1  =— f(  CjA'p+  2€PC|  )^6p 


(27) 


And  inserting  Equation(27)  into  Equation(ll)  and  simplifying  obtains: 


f(  C2A'p+  26PC,  )fcp=(l-f)<7y 


(28) 
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If  we  find  the  value  of  €p  from  Equation(28). 


i 


\ 


> 

\ 


i 

i 


€p  — _  YU,(a  _  a)A^  +  2^(1— f)<Ty  (  27  ) 

(  +:  5ep<0  -  Cooling,  — :  ^p>0  -  Heating  ) 

T 

If  this  value  of  €p  is  inserted  into  Equation(24)  and  (25),  the  value  of  €  and  the 
value  of  6-0  can  be  found,  and  if  these  values  are  inserted  into  Equaiion(5)  the  value 

_I  ij 

of  a  can  be  determined.  Therefore,  inserting  this  value  into  Equation(16)  we  are 
>j 

ultimately  able  to  obtain  total  average  stress  in  the  matrix.  Refer  to  Appendix(  0  ) 
for  the  constants  and  actual  values. 

3.  Creep 

Up  to  this  point  we  have  been  determining  the  internal  stress  both  in  the 
case  of  elastic  matrix  and  elastic  inclusion  and  in  the  case  of  elastoplastic  matrix 
and  elastic  inclusion.  But  in  these  two  cases  the  actual  effect  on  creep  deformation  is 
the  internal  stress  in  the  case  of  elastoplastic  matrix  and  elastic  inclusion. 
Accordingly,  by  inserting  the  value  of  the  internal  stress  into  Equation  (13),  we  can 
plot  the  microcreep  deformation  phenomenon  as  a  function  of  time  as  shown  in  Fig. 
3.2.  Refer  to  Chapter  IV  part  A  for  more  detailed  information. 


Ji  i 


IV.  RESULTS  AND  CONCLUSIONS 


A.  RESULTS 

The  thermo— mechanical  data  of  the  matrix  and  fiber  for  the  theoretical 

calculations  are  obrtained  from  the  [ref  13]. 

Annealed  2024  A1  matrix: 

Em  =  47.5  GPa 

a  =  47.5  MPa 

y 

v  =  0.3 

a=23.6*10_6/K  (31) 

SiC  fiber: 

Ei  =  427  GPa 
f  =0.2 

a  =  4.3X10-6/  K 

l/d=  1.5  (32) 

Where  the  average  value  of  the  fiber  aspect  ratio(l/d)  was  used  [ref  13].  The 
temperature  drop  aT  is  define  as 

aT  =  T1—  To  (33) 

Where  T1  is  taken  as  the  temperature  below  which  dislocation  generation  is 
minimal  during  the  cooling  process  and  To  is  the  room  temperature.  Thus,  for  the 
present  composite  system  aT  is  set  equal  to  — 200K.  From  the  data  given  by 
Equation  (31),  (32)  and  the  use  of  <  a >B  =  — we  have  computed  the  stresses. 
Next,  the  thermal  residual  stresses,  averaged  in  the  matrix  of  SiC  fiber/2024  AI,  are 


►  it* till 


r«*t  A'«»r »4  . 


predicted  by  (17)  and  the  result  on  <  r  >B  as  plotted  in  Fig.6,  where  33  denote 
the  component  along  the  longitudinal  direction(z),  The  average  theoretical  thermal 
residual  stress  is  predicted  to  be  tensile  in  nature,  and  the  average  residual  stress  in 
the  longitudinal  direction  to  be  larger  than  the  average  residual  stress  in  the 
transverse  direction  [Table  5].  The  fiber  aspect  ratio(of=l/d)  of  SiC  fibers  has  been 
observed  to  be  variable  [Appendix  A,B).  In  the  present  model  we  have  used  the 
value  of  1/d,  1.5  to  predict  the  thermal  residual  stress  of  the  composite. 


TABLE  3:  The  value  of  €~^  (  o=l,  1.5) 


TABLE  4:  The  value  of  6^  (a=l,  1.5) 


of  1/d) 

r1 

e  1  ■ 

T 

€ 

1 1 

2  2 

3  3 

1 

0.0076 

0.0076 

0.0076 

1.5 

0.0025 

0.0025 

0.0115 

; 

! 


TABLE  5:  The  value  of  (<*=  1,  1.5) 

»j 


a(  1/d) 

(r~\  MPa) 

fT^MPa) 

o'1  (MPa) 

1 1 

22 

33 

1 

-276.23 

-276.  2  3 

1.5 

-154.82 

-154.8  2 

TABLE  6:  The  value  of  <cr  ( cr=l,  1.5) 


o(l/d) 

<<7-i>(MPa) 

1 1 

<<r-1>(  MPa) 

22 

<tr_1>(MPa) 

33 

1 

55.245 

55.245 

55.24  5 

1.5 

30.905 

30.905 

70.998 

TABLE  7: 


a(l/d) 


Comparison  value  of  < 


(t~1  >a  with  TODAY'S 
'i 


OUR  MODELS  1  S 
fr~3  >  n  (  MPaj 


TAYA ' 

“ - =1“ 

<  ^33 


s _ 

>n(MPa) 


< 


Also  in  Fig.  4.1  we  observe  that  when  the  aspect  ratio(o)  increases,  the 
internal  stress  also  increases  with  it. 


Fig.  4.1  Internal  stress  vs  aspect  ratio 
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In  Fig.  4.3,  in  the  case  of  an  ellipsoid  (a=1.5),  the  values  of  internal  stresses 
increase  with  an  increase  in  volume  fraction. 


Fig.  4.3  Internal  stress  vs  volume  fraction(-f) 


Next,  if  we  compare  the  theoretically  predicted  value  of  the  internal  stress  with  the 
value  determined  by  R.  J.  ARSENAULT  and  M.  TAYA,  we  see  that  our  model's 
value  is  larger  than  the  value  which  is  obtained  by  R.  J.  ARSENAULT  and  M, 
TAYA,  as  shown  in  Table  7.  They  used  the  material  properties  as  follows: 

Annealed  6061  A1  matrix: 

E=47.5  GPa 
<r  =47.5  MPa 

y 

r/  =0.33 

a  =23.6*10-6/K 

SiC  whisker: 

E,=427GPa 
f  =0.17 

a*=4.3*10~6/K 

l/d=l.S 

Finally,  Fig.  3.2  is  a  graph  showing  the  strain  as  a  function  of  time,  when  the  aspect 
ratio  is  1.5  and  the  value  of  the  thermal  residual  stress  is  inserted  ir.to 
Equation!  13).  If  we  analyze  Fig.  2.3  we  can  see  that  the  creep  deformation 
phenomenon  is  due  to  the  internal  thermal  residual  stress  of  the  MMCs. 

B.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  object  of  this  research  is  to  obtain  the  value  of  the  thermal  residual  stress 
of  an  Al/SiC  composite  using  Eshelby's  theoretical  model,  and  then  to  determine 
what  effect  this  thermal  residual  stress  value  produces  on  creep  deformation. 
Because  exact  creeping  behavior  of  this  material  is  difficult  to  determine  when 
analyzing  creep  deformation,  Andrade's  model  was  used  with  properly  selected 


constant  values.  But  through  this  theoretical  approach  the  following  conclusions  can 
be  obtained. 


1.  Thermal  residual  stress  due  to  the  difference  in  thermal  expansion 

coefficients  may  be  estimated. 

a.  By  means  of  the  theoretical  model,  we  can  see  that  the  thermal 

residual  stresses  increase  when  the  volume  fraction  of  the  inclusions 
increases. 

b.  In  the  case  of  sphere  inclusions,  that  is,  when  the  aspect  ratio  o=l,  the 

lateral  stress  is  equal  to  the  longitudinal  stress. 

c.  In  the  case  of  an  ellipsoid  inclusions  where  cr=1.5  the  longitudinal  stress 

is  greater  than  the  lateral  stress. 

d.  The  thermal  residual  stresses  of  Ai/SiC  composites  increased  when 

the  value  of  the  aspect  ratio  of  the  inclusions  increases. 

2.  Microcreep  deformation  can  be  estimated  in  the  model  by  using  the 

thermal  residual  stresses  and  the  Andrade's  model  of  creep 
deformation 

a.  Microcreep  deformation  is  due  to  the  internal  thermal  residual 

stress  of  the  MMCs. 

b.  Dimensional  stability  of  components  will  be  influenced  by  the 

behavior  of  the  composite. 

3  Recommendations 

a.  Presently,  this  thesis  has  only  dealt  with  average  thermal  residual 

stress  from  a  overall  point  of  view,  and  we  need  more  detailed 
local  residual  stresses  surrounding  SiC  particles  should  be 
estimated  for  analyzing  creeping  behavior. 

b.  In  the  future,  research  should  be  done  concerning  relaxation  during 

cooling. 

c.  The  creep  behavior  used  in  the  current  model  is  the  Andrade's 

approx  -  .ation.  The  real  creep  deformation  for  this  Al/SiC 
composj.e  should  be  further  studied  in  order  to  get  better 
understanding  of  the  microcreep  of  the  material  based  on 
approximation  values. 


APPENDIX  A 


PROGRAM  FOR  VALUE  OF  Sjjkl(a=1.5  -  5) 
1.  PROGRAM(FORTAN) 


REAL*  8  NU ,  ALP ,  G ,  TMP  ,  Sllll ,  S3333 , 31122  ,  S1133 , S3311 
NU-0.3 

WRITE (6, 90)  i 

WRITE (6, 95) 

ALP-1.5 

300  TMP-ALP**2-1. 

C- (ALP* (ALP*DSQRT (TMP) -DLOG (ALP+DSQRT (TMP) ) ) ) / (TMP** (3 . /2 . ) ) 
Sllll- (3. *ALP**2) / (8 . • (1 .-HU) *TMP) 

♦■*• ( 1 . -2 . *NU-9 . / ( 4 . *TMP ) ) *G/ ( 4 . * ( 1-NU) ) 

S3333- ( (1.-2. *NU+ ( 3 . *ALP*  *2-1 . ) /TMP- (1.-2. *NU+3 . *ALP*  *  2/TMP ) 
+*C) ) / ( 2 . * (l.-NU) )  T  r 

S1122- (ALP **2/ (2. *TMP) - (1. -2*KU) -3*C/ (4 . *TMP) ) / (4 . * (l.-NU) ) 
S1133— ALP* * 2/  ( 2 .  *  ( 1 .  -HU)  *TMP )  ♦  ( 3 .  *ALP* *2/TMP-  ( 1 . -2 .  *NU) )  *G 
+/ (4. * (l.-NU) )  j 

S3311—  (1  .-2.  *NU+1 ./  (TMP) )  /  (2. •  (l.-NU) )  +  (1.-2.  *NU+3 ./  (2 .  *TMP)  ) 
+*G/(2.* (l.-NU) ) '  ! 

WRITE (6 . 100)  ALP, Sllll , S3333 , S1122 , S1133 , S3311 
ALP-ALP-*-.  1  I 

IP (ALP. GT. 5.)  GOTO  400 
GOTO  300 


90  FORMAT ( IX , '  ALP 

+  S3311 ' ) 

95  FORMAT (IX,'  — — 

+-«—') 

100  FORMAT (IX, 6 (F10. 5) ) 
400  STOP 
END 


Sllll 


S3333 


S1122  S1133 


2.  Sj.kl  VALUE  WITH  ASPECT  RATIO 


ALP 

sun 

S3333 

S1122 

S1133. 

S3311 

1.50000 

.58078 

.39378 

.01421 

.08396 

.01945 

1.60000 

.58828 

.37371 

.01603 

.08992 

.01623 

1.70000 

.59495 

.35518 

.01764 

.09549 

.01350 

1.80000 

.60089 

.33804 

.01906 

.10071 

.01117 

1.90000 

.60620 

.32215 

.02032 

.10560 

.00918 

2.00000 

.61098 

.30739 

.02145 

.11018 

.00747 

2.10000 

.61529 

.29366 

.02246 

.11448 

.00600 

2.20000 

.61920 

.28087 

.02337 

.11851 

.00472 

2.30000 

.62274 

.26893 

.02419 

.12231 

.00362 

2.40000 

.62598 

.25777 

.02493 

.12588 

.00266 

2.50000 

.62893 

.24732 

.02560 

.12924 

.00183 

2.60000 

.63164 

.23752 

.02622 

.13241 

.00111 

2.70000 

.63412 

.22832 

.02678 

.13540 

.00048 

2.80000 

.63641 

.21967 

.02729 

.13822 

-.00007 

2.90000 

.6385? 

.21152 

.02776 

.14090 

-.00056 

3.00000 

.64048 

.20384 

.02819 

.14342 

-.00098 

3.10000 

.64229 

.19659 

.02859 

.14582 

-.00134 

3.20000 

.64397 

.18974  i 

.02896 

.14809 

-.00167 

3.30000 

.64554 

.18326 

.02930 

.15025 

-.00195 

3.40000  • 

.64700 

.17713 

.02962 

.15230 

-.00219 

3.50000 

.64836 

.17131 

.02991 

.15425 

-.00241 

3.60000 

.64963 

.16579 

.03019 

.15610 

-.00259 

3.70000 

.65082 

.16054 

.03044 

.15787 

-.00276 

3.80000 

.65194 

.15555 

.03068 

.15955 

-.00290 

3.90000 

.65299 

.15081 

.03090 

v  .16115 

-.00302 

4.00000 

.65398 

.14629 

.03111 

.16269 

-.00312 

4.10000 

.65491 

.14197 

.03130 

.16415 

-.00321 

4.20000 

.65579 

.13786 

.03149 

.16555 

-.00329 

4.30000 

.65662 

.13393 

.03166 

.16689 

-.00335 

4.40000 

.65740 

.13018 

.03182 

.16818 

-.00340 

4.50000 

.65814 

.12659 

.03197 

.16940 

-.00345 

4.60000 

.65884 

.12316 

.03212 

.17058 

-.00348 

4.70000 

.65951 

.11987 

.03225 

.17171 

-.00351 

4.80000 

.66014 

.11671 

.03238 

.17280 

-.00353 

!4. 90900 

.66074 

.11369 

.03250 

.17384 

-.00355 

APPENDIX  B 


PROGRAM  FOR  VALUE  OF  Sijk](o=1.0,  «) 
1.  PROGRAM(MATLAB) 


* 

ASPHERE  ( ALPHA-1. 0) 

*  ! 

for  1-1:6; 

nu-i/10.;  j 

den-15.1*  (1-nu) ; 
ell (i) -<7-5*nu) /den; 
sl2(i)-(5*nu-l)/den; 

■13 (i) -sl2 (i) ; 
s31 (i) -sl2 (i) ; 

■33 (i)-sll(i) ;  .  • 

*  * 

VALPHA- INFINITY  » 

* 

den2-8* (1-nu) ; 

sail  ( i )  «•  (5-4.  *nu)/den2; 

8812 (i) -- (l-4*nu) /den2 ; 

ssl3 (i) -nu/ (2. * (1-nu) ) ; 

ss31(l)-0; 

sa33 (i)-0; 

end; 

re8ultl«[sll;sl2;8l3;831;s33;] ' 
result2- (sail ;ssl2; sa!3; ss31 ;ss33] ' 
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r.u 

r 

\ 

Sll 

S12 

S13 

S31 

S33 

resultl  - 

i 

C.l 

0.4815 

-0.0370 

-0.0370 

-0.0370 

0.4815 

0.2 

0.5000 

0 

0 

0 

0.5000 

0.3 

0.5238 

0.0476 

0.0476 

0.0476 

0.5238 

C.  4 

0.5556 

0.1111 

0.1111 

0.1111 

0.5556 

0.5 

0.6000 

0.2000 

0.2000 

0.2000 

0.6000 

C.5 

0.6667 

0.3333 

0.3333 

0.3333 

0.6667 

1 

result2  * 

1 

0.1 

0.6389 

-0.0833 

0.0556 

0 

0 

0.2 

0.6563 

-0.0312 

0.1250 

0 

0 

0.3 

0.6786 

0.0357 

0.2143 

0 

0 

0.4 

0.7083 

0.1250 

0.3333 

0 

0 

0.5 

0.7500 

0.2500 

0.5000 

0 

0 

0.6 

1 

0.8126 

0.4375 

0.7500 

0 

0 

APPENDIX  C 


PROGRAM  FOR  VALUE  OF  THERMAL  RESIDUAL  STRESS, 

p  T*  -C  -I  -I 

€j  .  ,  €  j  j  ,  6:  ~  .  (T  ,  <  (T  > 

1  1  1  ij  ij  n 

1.  PROGRAM(MATLAB) 

.r:q 

VThis  program  computes  stress  values  and  constants 
% 

%  Initial  values 

•  ^  — ■ 

nu-0 .  3 ; 

em-47.5e9; 

ei-427.0e9; 

sy-47.5e6; 

f-0.2; 

dt--20Q 

ap-23.6e-6; 

as-4.3e-6; 

mu-0 . ;  ’ r 

% 

%  Eshelby ' s  tensor  components 

% 

sll-0.6786; 

833-0.0000; 

812- 0.0357 ; 

813- 0.2143; 

&31-0.0000; 

%  . 

%  Constants  , 

\  ' 

app-(as-ap) *dt;  • 

g-em/ ( (l.+nu) * (1.-2. *nu) ) ; 

h-ei/(  (l.+nu)  •  (1.-2. *nu) )  ; 

j-(em-ei) /( (l.+nu) * (1.-2. *nu) ) ; 

a-f+sll* (l.-f ) ; * 

b-sl2* (l.-f) ; 

c-s!3* (l.-f) ; 

e-f+s33* (l.-f) ; 

d-s‘31*  (l.-f) ; 
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continued 


tl»rl+r2«’r3-ql+q2+q3; 

t2« (rl+r2-2*r3 )  /2 .  - (ql-q2+2*q3) /2 . ; 

t3*r4-*-r5*r6-p4-p5-p6; 

t4«(r4+r5-2*r6) /2.- (p4+p5-2*p6) /2. ;  * 

t5«r7+r8+r9-pl-p2-p3; 

t6- (r7+r8-2*r9)  /2.-  (pl+p2-2*p3)  /2.  ; 

xl-g«  ( (l.-nu)  *tl+nu* (t3+t5)  > ; 

x2»g* ( (l.-nu) *t2+nu« (t4+t6) ) ; 

x3»g*  ( (l.-nu)  •t3+nu« (tl+t5) ) ; 

x4»g*  ( (l.-nu)  *t4+nu*  (t2+t6) ) ; 

x5»g* ( (l.-nu) *t5+nu* (tl+t3) ) ; 

x6-g* ( (l.-nu) *t6+nu* (t2+t4) ) ; 

cl*( (x2+x4) /2.-x6) * (2/2. ) ; 

c2«( (xl*x3)/2.+x2+x4-x5+x6) *  It/2,) ; 

\ 

\  Temprature  dependent- parameter 

*  *  * 
ep— c2*app/(2.*cl)+(l.-f ) •sy/(2.*cl) ; 

% 

%  Elastic  moduli  of  the  matricies 

* 

aaa-g* [ (1-nu)  nu  nu  0  0  0 ;nu  (1-nu)  nu  0  0  0;nu  nu  (1-nu)  000; 
000  2.*nu  0  0;0  0  0  0  2. *mu  0;0  0  0  0  0  2.*mu]‘; 

t 

4 

Plastic  strain  'in  the  matrix 
ijp«[-ep/2  -ep/2  ep  0  0  0]' 

Uniform  transformation  strain  of  inclusion 
eijts"[app+ep/2  app+ep/2  app-ep  000]' 


rva  •  '*/Sgg33«  >&&&%• 


K.-T  MJTKX 


- . ■-  JJ '.wwm U w WTTfW U WWW W  L'W  JV W WIW1  ww IPJ'A  un^iunm wnnww 


continued 


eijt-tql* (app+ep/2) -q2* (app+ep/2)-q3* (app-ep)  . . 
p4* (app+ep/2) +p5* (app+ep/2) +p6* (app-ep)  . . 
pi*  (app+ep/2)+p2*  (app+ep/2) +p3*  (app-ep)  0  0  0] 
eij_c»[rl*  (app*ep/2)+r2*  (app+ep/2)+r3*  (app-ep)  . . 
r4*  (app+ep/2 )  +r5*  (app+ep/2) ♦rS*  (app-ep)  . . 
r7* (app+ep/2) +r8* (app+ep/2)+r9* (app-ep)  000) 

% 

%  Average  internal  stress  in  inclusion 

% 

si j_i®aaa* (eij_c-ei jt) 

% 

*  Average  stress  in  the  aatrix 

% 

si j»*-f *sij_i 
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2.  CONSTANTS  OF  ELASTIC  MATRIX  AND  ELASTIC  INCLUSIONS 


B=Sn(  1—1 )  +  i 
D=S,  3(l-f) 

F=S3.(l-f) 

H=Ei/(l+i/)(l-2i^) 

K= J((  l-t/)B + i/C+ i/F  H  l-i/)G 
M = J(  ( 1— i')D + i'D + iTS )— j'G 

PsJ((MB+iCt«'FHMG 
R ..= J( ( ( 1  —  v )  F + vB + i'C  )—i'G 
V=J((l-i/)B+2i/DHl-i/)G 
V=(QK  -  NM)/K 
X=(SK  -  RL)/K) 
Z=a'h(1+z')(R-K)/K 

v'=(w-x'v)/u 


L— Si2l  1—1 )  +  1 
E=S33(1— f)  +  1 
G  =  Em/(l+z/)(  l-2i/) 

J=(Em— Ei)/(l+z^)(l— 2z/) 

L= J((  1— i/)C+i/B+t/F)— zX3 
N =J(  ( 1— i/)  D + i/D + i'E  )—i'G 
Q=.J(fl-,/)D+i/D+i/E)-i/G 
S=  J((l— i')F+z'C+i'B)— i'C 
U=(PK  —  NL)/K 
W=a'h(1+j/)(N-K)/K 
Y=(TK  -  MR)/K 
X ' =(  ZU-X  W )/( U  Y-X  V ) 

z^i-y'l-x'm-ha'ii+^d/k 


PI! OGHAM  FOR  CREEP  REFORMATION 


PIlOGRAM(  MATLAB) 


ep0*O . 001 ; 
became . 0465 ; 
k *7.2e-7; 
for  tt*l : 91 ; 
kv ( tt) *tt; 

t=1.0«10~ ( ( itt-l)/30.) )  ; 

eep(tt)=ep0«  ( 1 . fbeta" t ~  (1/3)  )  «exp(k*t)  ; 

end; 

plot (kv, eep) 
semilogx; 

xlabel ( ‘Time  t *), ylabel (' Strain  eps * ) 
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